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ABSTRACT 


This  report  describes  the  results  of  an  empirical  study 
to  develop  a  normalizing  transformation  for  sample  multiple 
coherence.  The  "Fisher  z-transformation"  is  employed.  The 
expected  value  (including  a  bias  term)  and  variance  of  the  trans¬ 
formation  have  been  experimentally  determined.  Numerical 
^values  of  the  transformation  which  is  developed  (including  the 
bias  term  and  variance)  may  be  obtained  with  a  reasonable 
amount  of  computation.  Tables  of  the  Gaussian  distribution 
can  then  be  used  to  obtain  confidence  limits  and  perform  sta¬ 
tistical  tests. 

The  computational  methods  and  the  digital  computer 
program  used  for  the  study  are  described  in  detail.  Flow 
charts  of  the  program  are  given.  Numerical  results  from 
the  program  results  are  presented.  Examples  of  the  use  of 
the  transformation  are  given  for  developing  confidence  limits 
for  multiple  coherence.  A  recommendation  for  a  further 
theoretical  study  is  presented. 
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1.  INTRODUCTION 


In  measuring  coherence  when  moderately  long  or  long  records  are 
available,  one  often  has  the  "effective  degrees -of-freedom"  parameter  n  >  20. 

A  large  set  of  tables  consisting  of  nine  voliunes  has  been  developed  by  the 
Rocketdyne  Division  of  North  American  Aviation.  These  tables  provide  the 
sampling  distribution  of  coherence  when  the  effective  degrees -of-freedom  is 
roughly  less  than  20,  and  when  the  ntimber  of  variables  involved  in  the 
coherence  relation  is  less  than  eleven.  This  set  of  tables  is  listed  as 
Reference  1  in  this  report. 

In  regard  to  many  practical  applications,  it  is  desirable  to  have  a  tract¬ 
able  and  reasonably  accurate  approximation  to  the  distribution  of  sample 
coherence  for  n  >  20.  Such  an  approximation  is  especially  needed  for  the  range 
20  <  n  <  150.  The  degrees -of-freedom  parameter  in  practical  spectral  analysis 
often  falls  in  this  range,  and  hence  the  degrees -of-freedom  for  coherence 
functions  obtained  from  the  spectral  analysis  also  falls  in  this  range. 

The  object  of  the  study  described  in  this  report  was  that  of  obtaining  a 
Gaussian  distributed  approximation  to  sample  coherence.  A  reasonable  approxi¬ 
mation  which  has  a  roughly  Gaussian  distribution  has  been  empirically  obtained  by 
a  digital  computer  numerical  analysis.  This  approximation  appears  to  be  useful 
in  the  range  where  true  coherence  varies  from  roughly  0.4  to  0.95,  and  the 
degrees- Qf-freedom  parameter  is  larger  than  20.  A  slightly  more  complicated 
version  of  the  transformation  involved  in  the  approximation  gives  corrections 
which  extend  the  usefulness  of  the  approximation  down  to  approximately  a  true 
coherence  value  of  0.1,  and  to  the  degrees -of-freedom  parameter  as  low  as  10. 

The  approximation  developed  is  very  useful  in  that  the  need  to  have  a  large  set 
of  tables  available  for  many  practical  applications  is  eliminated.  One  instead 
can  substitute  a  table  of  the  Gaussian  distribution  function  which  is  widely  and 
conveniently  available.  A  certain  amount  of  computation  is  necessary  in  order 
to  transform  the  value  of  sample  coherence  to  the  approximately  Gaussian  distri¬ 
bution  variable:  however,  the  amount  is  not  unreasonable.  The  transformation  and 
the  empirical  procedures  used  in  obtaining  the  final  version  will  now  be  described. 

1 


2.  THE  NORMALIZING  TRANSFORMATION 


The  probability  density  function  for  sample  coherence  is  given  by  the 
following  general  formula: 


C(x|  n,  p, 


Y^)  = 


r(n) 


r(p-l)r(n-p+l) 


2  n  d-2  n-D  2 

(1  -  Y  )  (1  -x)  ^F(n,  n;p-l;Y  x)  (1) 


(0  <  X  <  1) 


where 


n  = 


P  = 
2 


X  = 


Y  = 


F(n,  n;p-l  ;y  x)  = 


effective  number  of  degrees-of-freedom  (d.  f. ) 

effective  number  of  records 

true  coherence 

sample  coherence 

hypergeometric  function 


A  closed  form  expression  for  the  cumulative  distribution  function  corresponding 
to  Eq.  (1)  is  known.  The  closed  form  expression  is,  however,  intractable  for 
n  in  the  range  desired.  That  is,  the  computational  requirements  to  obtain  a 
numerical  value  for  a  specific  n  requires  more  than  a  reasonable  amount  of 
computation.  Recursion  relations  on  n  have  also  been  developed.  The  tables 
described  in  Reference  1  were  computed  using  these  relations.  These  recursion 
formulas  are  suitable  for  obtaining  results  for  n  +  1  in  terms  of  the  results  for 
n.  However,  they  are  not  suitable  for  directly  computing  a  numerical  value  for 
any  arbitrary  value  of  n. 

The  asymptotic  (n  approaching  infinity)  distribution  corresponding  to 
Eq.  (1)  has  been  investigated.  However,  comparisons  with  the  existing  tables 
of  Reference  1  indicated  that  approximations  based  on  such  asymptotic  theory 
would  be  far  from  reasonably  accurate  for  20  <  n  ^150.  Other  preliminary 
studies  have  indicated  that  simplicity  and  reasonable  accuracy  may  be  achieved 
by  a  "normalising”  (Gaussian  producing)  transformation  of  the  form  described 
below. 
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Let 


z=-in 


1  tYI 

r  yir 


1=  tanh  ^ 


(2) 


tanh  ^  (\) 


(3) 


In  Eqs.  (2)  and  (3),  x  =  ^  =  'Y  sample  coherence  and  y  =  true  coherence. 

The  random  variable  z  is,  to  a  close  approximation,  normally  (Gaussian) 
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distributed  N(p,  cr  ).  The  mean  p  is  given  by  an  expression  of  the  form 


iLi 


n  -  p  +  k 


=  C  +  b(Y) 


(4) 


and  the  variance  by 


2 

<r  = 


1 


2(n  -  p  +  k) 


(5) 


In  Eqs.  (4)  and  (5),  the  constant  k  and  the  function  f(y)  are  not  yet  specified. 

The  choice  of  the  form  of  the  transformation  given  by  £q.  (2)  is  motivated 
by  a  comparison  of  multiple  coherence  with  the  multiple  correlation  coefficient 
of  classical  statistics.  When  this  transformation  (termed  the  'Fisher  z> 
transformation")  is  applied  to  the  ordinary  correlation  coefficient,  one  obtains 
a  variable  which  is  normally  distributed  to  a  very  high  degree  of  accuracy. 

Biases  do  exist  in  the  transformation  of  the  ordinary  correlation  coefficient,  but 
they  are  usually  negligible  for  practical  purposes.  However,  for  the  case  of 
multiple  coherence,  the  transformation  by  itself  is  not  quite  so  suitable.  The 
bias  term  denoted  by  b{y)  can  be  significant.  One  of  the  major  objects  of  the 
study  described  in  this  report  was  to  determine  the  numerical  behavior  of  this 
function  b(y)  and  to  determine  if  it  could  be  represented  by  a  reasonably  simple 
form. 
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The  asymptotic  (n  approaching  infinity)  distribution  of  the  random  variable 
z  is  normal  with  a  mean  of  t  and  a  vai  ance  of  l/2n,  provided  y  ^  0.  Initial 
comparisons  which  were  made  with  tabulated  results  indicated,  however,  that 
for  n  near  20,  the  asymptotic  distribution  is  still  a  far  from  acceptable 
approximation.  One  may  view  Eqs.  (4)  and  (5)  as  giving  finite  corrections  to 
the  asymptotic  distribution  of  the  random  variable  z.  For  various  hand  calcu¬ 
lated  examples,  corrections  of  the  form  indicated  by  Eqs.  (4)  and  (5)  were  intro¬ 
duced  and  acceptable  approximations  to  the  distribution  of  sample  coherence  were 
obtained.  These  hand  calculations  motivated  the  performance  of  a  more  compre¬ 
hensive  examination  by  means  of  a  digital  computer  program. 

The  object  of  this  program  was  to  determine  if  z  is  a  useful  transformation 
tor  coherence  when  suitable  bias  and  variance  terms  were  introduced.  Thus, 
determining  the  form  of  the  bias  b(y)  and  the  value  of  k  which  would  fit  in  the 
expressions  for  the  mean  and  variance  were  to  be  determined  by  the  digital 
computer  program.  The  next  section  of  the  report  describes  the  computer 
program  which  was  written,  and  the  section  after  that  presents  the  results  which 
were  obtained. 
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3.  COMPUTATIONAL  PROCEDUKES 


i 


The  basic  method  of  determining  the  accuracy  of  the  normalizing 
transformation  was  to  compute  results  using  the  transformation  and  compare 
them  with  the  existing  tables  of  sample  coherence.  These  tables  are  avail* 
able  for  n  up  to  roughly  20.  It  is  assumed  that  if  an  acceptable  transformation 
is  determined  for  this  small  n,  then  the  transformation  will  also  be  accept¬ 
able  for  larger  n.  It  was  decided  to  determine  tlie  bias  of  the  transformation 
by  comparing  the  deviation  corresponding  to  the  50%  point  (median)  of  the 
sample  coherence  table  from  that  corresponding  to  the  true  value  of 
coherence.  It  was  also  decided  to  establish  the  variance  for  the  transfor¬ 
mation  by  comparing  the  95%  point  of  the  distribution.  That  is,  for  a 
Gaussian  distribution  the  95th  percentile  should  be  1. 645  standard  deviations 
away  from  the  mean.  Therefore,  the  difference  between  the  50th  and  95th 
percentiles  was  set  equal  to  1.645  standard  deviations.  This  allowed  the 
determination  of  the  variance. 

After  the  bias  and  the  variance  were  determined,  comparisons  for 
several  probability  points  were  then  made  of  the  hypothesized  normally 
distributed  transformed  values  with  the  exact  tabulated  distribution  values. 

To  be  more  specific,  the  computational  steps  proceed  as  follows. 

(1)  Select  a  value  of  true  coherence. 

(2)  Search  the  appropriate  table  of  sample  coherence  for  the  50th 
percentile  of  the  distribution  of  ^  given  y^.  This  procedure 
actually  requires  interpolation  in  the  existing  table.  This  inter¬ 
polation  was  performed  using  3  point  quadratic  interpolation 
defined  by  the  equation 


A^  tP-P^XP-P^) 
•SO  (Pi-P^XPj-Pj) 


4- 


(p-pjHp-Pj) 

(P^-PjHp^-Pj) 


4 


(p-PjHP'P^) 

(pj-p,«pj-pj) 


(6) 


5 


JP 


In  Eq.  (6), 


p  =  .50 

Pi  I  Pti  Po  three  distribution  function  values 

X  M 

which  bracket  p  =  ,  50 

arc  the  arguments  corresponding  to  tlie 
distribution  function  values  p^  ,  p^  f  and  p,- 

The  bias  at  the  point  is  then  determined  by  transforming 
th'=  50%  point  of  sample  coherence  and  comparing  that  with  the 
transformed  value  of  the  true  coherence.  That  is,  b  is 
defined  by  the  equation: 

^0.50  =  ^0.50  ’ 


=  Vso  -  ^ 


After  the  bias  is  determined,  then  the  variance  is  determined 

by  a  similar  procedure.  First,  the  95th  percentile,  y  nc  » 

.95 

obtained  by  interpolation  with  a  formula  the  same  as  Eq.  (6), 
but  with  the  following  quantities  replacing  those  in  Eq.  (6). 


p  =  .  95 

P/I  •  Pc  »  P/^  three  distribution  function  values  which 

4  b  D 

bracket  p=.  95  (these  replace  p. ,  i=l,2,  3  in 


Eq.  (6). 


^4  ’  ^5  ’  ^6  arguments  corresponding  to  p^  ,  p^ ,  p^ 


and  replace  »  i  =  1,  2,  3. 


6 


The  variance  for  the  transformed  variable  is  then  determined 
from  the  following  relation 


(8) 


Paints  on  the  normal  distribution  function  with  a  mean 


fi.  =  t  +  b  =  z  (9) 

'  0.50 

and  a  variance  determined  by  Eq.  (8)  are  next  computed  for 
values  of  in  the  range  0.  05  to  0.  95  in  steps  of  0.  05. 

(6)  Since  the  asyniptotic  form  of  the  variance  of  the  sampling 
distribution  of  the  transformed  coherence  is  of  the  form 
0-^  =  l/2n,  it  was  decided  to  obtain  an  ’’effective"  number  of 
degrees -of-freedom  n  from  the  relation 


n 

e 


(10) 


(7)  Finally,  the  probability  values  which  have  been  computed,  the 
bias,  and  the  effective  degrees-of-freedom  are  printed  out. 


For  a  given  degrees-of-freedom  parameter  n  and  a  given  number  of 

variables  p,  the  above  procedure  was  repeated  for  values  of  true  coherence 

ranging  from  0.  00  to  0.  90  in  steps  of  0.  10,  and  from  0.  91  to  0.  99  in  steps 

of  0.  01.  It  was  also  suspected  that  the  results  would  vary  roughly  as  a 

function  of  the  square  root  of  n  .  Hence,  24  cases  were  run  with  V  n 

e  e 

running  from  roughly  2.  0  to  4.  5  in  steps  of  about  0.  5.  The  cases  w'  ich 
were  actually  run  are  given  in  Table  1  below.  In  Table  1,  n^  =  n-p  +  1. 
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TABLE  1 

Values  of  n  and  p  for  Which  Computer  Computations  Were  Performed 


1.  7 

2.4 

3.  0 

3.  5 

4.0 

4.4,  4.  4,  4.  3 

n 

e 

3 

6 

9 

12 

16 

19,  19,  18 

n,  p 

m 

B 

10,  2 

13,  2 

17,  2 

— 

5.  3 

8,  3 

11.  3 

14,  3 

18,  3 

21.  3 

■■■ 

9.  4 

12.  4 

15,  4 

19,  4 

22,  4 

n 

10,  5 

13,  5 

16,  5 

20,  5 

22,  5 

This  selection  of  values  of  n  and  p  provide  a  coverage  of  the  existing  tables 
for  2  <  p  <  5  and  allowed  for  the  investigation  of  the  behavior  of  the 
approximation  under  a  reasonably  wide  range  of  values  of  n. 

A  flow  chart  which  outlines  the  computational  procedure  used  is  given 
in  Figures  1  and  2.  The  polynomial  approximations  given  in  these  figures 
for  the  Gaussian  distribution  function  were  obtained  from  Reference  2.  These 
approximations  were  checked  with  existing  table  values,  and  agreement  was 
found  out  to  at  least  five  decimal  places  in  all  cases. 
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Figure  1.  Over-all  Computation  Procedure 
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Figure  2.  Interpolation  and  Gaussian  Distribution  Si'broutines 
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4.  NUMERICAL  RESULTS 


As  a  general  statement,  the  comparison  of  the  probabilities  obtained 
with  the  Gaussian  approximation  (transformation)  with  the  bias  and  variance 
adjustments  incorporated,  agreed  with  the  existing  multiple  coherence  tables 
to  two  or  three  places  for  fairly  large  degrees-of-freedom.  For  smaller 
degrees -of-freedom,  the  probabilities  were  still  often  quite  close  except  at 
times  toward  the  tails  of  the  distribution.  Even  here  the  errors  were  seldom 
worse  in  magnitude  than  0.  01  or  0.  02.  The  acbial  computer  output  of  all  the 
probability  values  computed  is  very  extensive,  and  the  complete  results  are 
not  included  in  this  report  in  the  interest  of  economy  of  space.  However,  the 
bias  values  obtained  for  all  the  cases  are  plotted  in  Figure  3  through  Figure  8 
(located  at  the  end  of  this  section)  along  with  two  approximations  to  the  bias 
correction  described  below.  Also,  an  example  of  the  agreement  in  probabilities 
is  given  later  in  this  section. 

An  analysis  of  the  numerical  values  obtained  for  the  effective  degrees- 
of-freedom  n^  obtained  via  the  variance  fitting  computations  indicates  that 
a  reasonable  value  to  use  for  n^  is  given  by 


and  the  variance  is 


n 


1 


n 

e 


p  +  2 


p  =  2 
p  >  2 


l/2n 

e 


(H) 


(12) 


This  value  becomes  less  accurate  as  coherence  and  degrees-of-freedom 
become  smaller  and  as  p  becomes  larger.  However,  it  is  a  quite  satis¬ 
factory  approximation  in  the  range  of  roughly  true  coherence  between  0.  4 
and  0.  95. 
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A  very  simple  straightforward  constant  bias  correction,  that  is,  one 
which  is  not  a  function  of  the  true  coherence  is  given  by  the  formula 


b  = 


P  -  ^ 

2(n  -  p  +  1) 


(13) 


This  value  is  plotted  in  Figures  3  through  8  as  the  constant  line.  From  this 
it  can  be  seen  that  for  fairly  large  degrees-of-freedom  this  value  of  b  is  a 
qi’ite  satisfactory  approximation  for  the  bias  for  true  coherence  values  of 
roughly  0.  4  through  0.  95. 


One  comment  should  be  made  about  the  plots  of  the  bias  correction 
at  this  time.  From  the  general  theory  of  the  approximations,  one  would 
expect  the  singularity  which  appears  to  occur  near  true  coherence  equal  to 
zero.  However,  there  is  no  immediate  explanation  for  the  apparent  singu¬ 
larity  which  occurs  at  true  coherence  equal  to  one.  This  leads  one  to 
believe  that  the  erratic  results  obtained  as  true  coherence  approaches 
closely  to  one  may  actually  be  due  to  computational  noise.  That  is, 
certain  computational  instabilities  could  possibly  be  occurring  in  this  area 
which  would  result  in  the  unexpected  results  obtained  here.  This  problem 
has  not  yet  been  resolved.  Further  theoretical  studies  will  tend  to  shed 
more  light  on  this  question.  A  modified  bias  correction  gives  somewhat 
better  results  for  lower  values  of  true  coherence  and  lower  degrees-of- 
freedom.  The  correction  is  given  by  the  formula 


"  2(n  *p+l)  ^  (Y  +  . 


1) 


P  •  ^ 


4(n-p  +  1) 


1  + 


yt.l 


(14) 


This  equation  is  plotted  as  the  dotted  lines  on  Figures  3  through  8.  Note 
that  this  bias  correction  has  the  disadvantage  of  having  true  coherence 
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appear  in  the  formula.  The  true  coherence  will,  in  general,  not  be  known 
and  hence  one  would  have  to  use  either  an  estimated  value  or  a  value  of  the 
sample  coherence.  The  use  of  this  formula  would  still  be  more  likely  to 
provide  a  slightly  better  bias  correction  term  than  that  given  by  £q.  (13). 
Again,  further  theoretical  studies  concerning  normalizing  transformations 
for  coherence  would  tend  to  resolve  this  question. 

The  plots  given  in  Figures  3  through  8  depict  empirically  the  behavior 
of  the  bias  term  h{y)  defined  by  Eq.  (4).  Two  approximating  forms  for  b(Y) 
have  been  given  by  Eqs.  (13)  and  (14).  It  is  this  empirical  behavior  whose 
analytical  form  should  be  determined.  The  computer  output  for  n  =  (n-p+k) 
as  defined  by  the  variance  formula,  Eq.  (5)  is  not  presented.  Equation  (11) 
summarizes  the  values  of  k  which  were  determined  as  suitable  approximations 
from  the  computer  study. 

A  typical  example  is  illustrated  in  Table  2  of  the  agreement  between 
the  sampling  distribution  probabilities  obtained  in  three  ways:  (a)  the  approxi¬ 
mation  using  the  bias  correction  of  Eq.  (13)  and  the  variance  of  Eq.  (12), 

(b)  the  digital  computer  fitted  approximation  using  the  median  bias  (see 
Figure  6)  and  the  variance  determined  from  the  95th  percentile,  and  (c)  the 
exact  probabilities  obtained  from  the  tables  of  Reference  1.  The  case  chosen 
is  for  =  0.  70.  n  =  15,  and  p  =  4. 
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TABLE  2 

Comparison  of  Approximate  and  Exact  Sampe  Coherence  Probabilities 


a2 

Y 

arc  tanhv 

Z 

(a)  P(2) 

(b)  P(z) 

(c)  P(z) 

.  50 

.  707 

.881 

-2.321 

.  010 

.008 

.  008 

.  55 

.  742 

.955 

-1. 944 

.  026 

.  022 

.  022 

.60 

.  775 

1.043 

-1.495 

.  068 

.  054 

.  054 

.65 

.  806 

1.117 

-1. 117 

.  132 

.  123 

.  123 

.  70 

.  837 

1.  211 

-  .  628 

.261 

.251 

.252 

.  75 

.866 

1.317 

-  .097 

.461 

.456 

.456 

.  80 

.  895 

1.447 

.  566 

.714 

.  710 

.707 

.  85 

.  923 

1.603 

1. 362 

.  913 

.916 

.910 

.90 

.950 

1.832 

2.  531 

.994 

.995 

.992 

=  0.  70,  n  =  15,  p  =  4 

(a)  Gaussian  Probabilities  Obtained  Using  Eqs.  (12)  and  (13) 
for  Variance  and  Bias.  Bias  =  b  =  3/24  =  .  125  , 

Standard  Deviation  =  a  -  =  .  196 

(b)  Gaussian  Probabilities  Obtained  Using  Digital  Computer 
Fitted  Median  Bias  and  95th  Percentile  Variance 

(c)  Exact  Probabilities  Obtained  Using  Tables  of  Reference  1 
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15 


Figure  3.  Graphs  of  Bias  Error  for  Normalizing  Transformation  for 
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Figure  4.  Graph*  of  Bia*  Error  for  Normalizing  Trantformation  for 


Constant  Lines 
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Figure  S.  Graphs  of  Bias  Error  for  Normalising  Transformation  for 


Figure  7.  Graphs  of  Bias  Error  for  Normalizing  Transformation  for 
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Figure  6.  Graphs  of  Bias  Error  for  Normalizing  Transformation  for 
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Figure  7.  Graphs  oi  Bias  Error  (or  Normalizing  Transformation  for 


Constant  Lines 
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Figure  8.  Graphs  of  Bias  Error  for  Normalizing  Transformation  for 


5.  CONFIDENCE  BOUND  EXAMPLES 


The  transformation  which  is  a  reasonably  accurate  normal  approxi- 

2 

mation  for  moderate  values  of  true  coherence  (0.4  <  y  <  0,95)  and  also 
useful  for  more  extreme  values  is  given  by  Eq.  (2). 

^  ,-lA  1,  1+Y 

z  =  tanh  Y  =  —  in  - 

2  ,  A 

1  -  Y 

where  ^  is  the  positive  square  root  of  sample  coherence.  In  Section  4, 
it  has  been  empirically  determined  that  z  is  approximately  normal  with  an 
approximate  mean  value 


1.  1+Y  P“1 

M-  =  T  ill  - - TTT 

2  1  -  Y  ”  P  +  1) 


where  y  is  ^be  square  root  of  true  coherence.  The  variance  of  z  is 
approximately 


^z  2(n-p  +  k)  ' 


k  =  1  for  p  =  2 
k  =  2  for  p  2 


From  the  above  equations,  if  one  selects  a  confidence  level  {I  -  a) 
and  is  given  a  value  of  y^  ,  then 


PjjL-aZ  ,  <z<|i  + 

■^z  z  l-cr/2  -  '^x 


*^z^l  -or/ 2  I 


1  -  a 


where  is  the  (1  -  a/ 2)  percentile  of  the  normal  distribution. 

Let  be  broken  into  two  parts 


t  =  yJn 

2  1  -  Y 


2(n  -  p  +  1) 
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Then,  Eq.  (17)  can  be  rearranged^  to  give  (1  -  a)  confidence  bounds  on  y 
which  are  defined  by  the  inequality 

tanh  (z  -  b  -  atZ^^  £  V  <  (z  -  b  + 

To  employ  the  foregoing  equations  in  developing  confidence  bounds  for 
the  various  types  of  coherence,  certain  adjustments  must  be  made  in  the 
parameters  n  and  p.  The  adjustments  are  discussed  in  Reference  3,  and 
will  be  reviewed  here.  First,  define  p'  =  the  total  number  of  variables  being 
considered,  n'  =  BT  as  the  number  of  degrees -of-freedom  in  the  associated 
spectral  analysis.  (In  many  documents,  the  number  of  degrees-of-frv’edom  is 
defined  as  n  =  2BT  exactly  double  the  value  used  here.  One  must  use  care  in 
comparing  degree-of-freedom  parameters  from  such  documents  with  this 
report.  ) 

(1)  For  general  multiple  coherence  between  one  output  variable 
and  p'  -  1  input  variables,  set  n  =  n'  and  p  =:  p'. 

(2)  For  marginal  multiple  coherence  between  one  output  variable 
and  p'  -  1  -  q  input  variables  (q  variables  being  neglected  in 
the  relation),  set  n  =  n*  and  p  =  p'  -  q. 

(3)  For  conditional  (partial)  multiple  coherence  between  one 
"conditioned"  output  variable  and  p'  -  1  -  i  "conditioned"  input 
variables  (the  linear  effect  of  i  variables  has  been  subtracted 
out),  set  n  =  n'  -  i  and  p  =  p'  -  i. 

(4)  For  marginal  conditional  multiple  coherence  between  one  con¬ 
ditioned  output  variable  and  p'  -  1  -  i  -  q  conditioned  input 
variables  (the  linear  effect  of  i  variables  has  been  subtracted 
out  and  q  variables  are  being  neglected),  set  n  =  n'  -  i  and 

p  =  p'  -  i  -  q. 

For  example,  for  ordinary  coherence  between  one  input  and  one 
output,  set  p  =  2  and  n  =  n'.  For  partial  coherence  between  one  input 
and  one  output  where  the  effects  of  a  second  input  are  subtracted  out,  set 
p  =  2  and  n  =  n'  -  1 . 
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Note  that  the  effect  of  the  adjustment  in  the  dewrees-of-freedom  will 

be  small  for  small  p'  and  large  for  n'  ,  say  n'  >  30.  The  bias  correction 

is  also  minor  for  many  practical  purposes  for  small  p'  and  large  n'.  For 

2 

example,  for  y  =  0.  9  [y  0.  95),  p  =  2,  and  n  =  25,  the  relative  bias  is 


b  0. 02 
t  '  1.47 


0.014 


1.4% 


The  bias  increases  at  least  linearly  with  p.  Hence,  for  larger  p,  the 
percentagewise  effect  becomes  significant. 

The  use  of  Eq.  (18)  is  illustrated  in  Figure  9  and  Table  3.  The 
99%  confidence  limits  for  n  =  BT  =  50  and  p  =  2  are  plotted  and  tabulated 
there.  The  dashed  lines  for  the  confidence  limits  where  is  less  than 
0.  3  is  to  indicate  that  the  approximation  is  becoming  unreliable  in  this  area. 
In  order  that  the  general  shape  of  the  approximate  and  true  confidence  limit 
curves  may  be  compared,  true  values  are  presented  in  Figure  10  for 
n  =  21  and  p  =  2.  The  values  for  these  curves  are  obtained  from  the 
tables  of  Reference  1. 
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Figure  9.  99%  Conlidence  Limits  for  Ordinary  and  Partial  Coherence; 

n  =  50,  p  =  2.  Obtained  by  Using  Equation  (18) 


TABLE  3 


Confidence  Limits  for  Ordinary  and  Partial  Coherence 


a2 

V 

.2 

.  3 

D 

.5 

D 

B 

.8 

H| 

liilDuiSjillilillli 

.  31 

.49 

.  58 

.66 

.  74 

.81 

.87 

.94 

IBSjS9HIh 

.  11 

.  19 

.29 

.41 

.54 

.67 

.83 
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Figure  10.  Exact  99%  Confidence  Limits  for  Ordinary  and  Partial  Coherence; 

0^21,  p  =  2  Stained  from  Tabulated  Oietribntlon  Function  of 
Reference  1 
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6.  CONCLUSIONS  AND  RECOMMENDATIONS 


A  useful  normalizing  transformation  for  sample  multiple  coherence 
has  been  empirically  developed.  The  computational  effort  required  to 
perform  the  transformation  is  reasonable  and  allows  the  use  of  Gaussian 
probability  distribution  tables. 

Several  aspects  of  the  transformation  are  not  completely  resolved 
by  this  empirical  approach.  The  transformed  variable  is  biased.  The 
numerical  behavior  of  this  bias  has  to  a  certain  extent  been  determined  as 
a  function  of  degrees  -of-freedom  n  and  number  of  variables  p.  This 
behavior  can  be  approximated  in  a  manner  which  is  satisfactory  for  multiple 
coherence  roughly  between  0.  40  and  0.  95,  However,  unexplained  numerical 
behavior  of  this  bias  occurs  when  multiple  coherence  is  close  to  1.  0.  Also, 
as  multiple  coherence  approaches  zero,  the  bias  seems  to  exhibit  singular 
behavior  which  is  more  or  less  expected  from  theoretical  considerations. 

The  variance  of  the  transformed  variable  was  empirically  developed 
as  a  result  of  this  numerical  study.  This  experimentally  determined  variance 
is  useful  for  about  the  same  ranges  of  the  parameter  values  as  for  the  bias. 

A  theoretical  study  to  attempt  to  determine  precise  analytical 
behavior  of  various  transformations  of  sample  multiple  coherence  would  be 
a  fruitful  undertaking.  A  study  of  this  sort  would  allow  more  confidence  to 
be  placed  in  the  reliability  of  transformed  values.  The  functional  forms  of 
the  bias  and  variance  could  be  determined  to  extend  the  usefulness  of  trans¬ 
formations  to  a  wider  range  of  parameter  values. 

The  multiple  coherence  function  is  closely  related  to  the  multiple 
correlation  coefficient  of  classical  statistics.  Hence,  it  is  expected  that  a 
detailed  theoretical  study  would  parallel  that  performed  by  Hotelling  for  the 
correlation  coefficient  (Reference  4).  There  seem  to  be  many  complications 
which  enter  into  an  examination  of  coherence  that  do  not  occur  for  the 
correlation  coefficient.  However,  the  results  of  such  a  study  for  coherence 
would  provide  useful  information. 
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